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Land topography can affect air radiation dose rates by locating radiation sources closer to, or 
further from, detector locations when compared to perfectly ffat terrain. Hills and slopes can also 
shield against the propagation of gamma rays. To understand the possible magnitude of topographic 
effects on air dose rates, this study presents calculations for ambient dose equivalent rates at a range 
of heights above the ground for varying land topographies. The geometries considered were angled 
ground at the intersection of two planar surfaces, which is a model for slopes neighboring flat land, 
and a simple conical geometry, representing settings from hilltops to valley bottoms. In each case 
the radiation source was radioactive cesium fallout, and the slope angle was varied systematically 
to determine the effect of topography on the air dose rate. Under the assumption of homogeneous 
fallout across the land surface, and for these geometries and detector locations, the dose rates at 
high altitudes are more strongly affected by the underlying land topography than those close to 
ground level. At a height of 300 m, uneven topographies can lead to a 50 % change in air dose 
rates compared to if the ground were uniformly flat. However, in practice the effect will more often 
than not be smaller than this, and heterogeneity in the source distribution is likely to be a more 
significant factor in determining local air dose rates. 


I. INTRODUCTION 

Surveys of gamma air radiation levels play an impor¬ 
tant role in quantifying radiological hazards if fallout 
from a nuclear accident or test contaminates an envi¬ 
ronment with radioactive material. The results typi¬ 
cally quoted from radiation surveys are values of the 
ambient dose equivalent rate, iL*(10) [pSvh“^], at 1 m 
above ground level, commonly referred to as the “air dose 
rate” [1]. 

Data obtained from air dose rate surveys can be used to 
determine the type and quantity of radioactivity present 
in the environment. The radionuclides in the fallout can 
be identified by analyzing the energy spectrum of the 
gamma radiation at the detector [2, 3]. The activity of 
radionuclides within the ground can then be estimated by 
employing conversion coefficients between 1 m ambient 
dose equivalent rates and soil activity levels [4]. 

With this knowledge, it is possible to calculate how 
air dose rates will change in future by performing decay 
calculations and modeling radionuclide migration [5, 6]. 
It is also possible to estimate external gamma radiation 
doses for exposed populations, and internal doses due to 
inhaling re-suspended activity or ingesting contaminated 
food-stuffs. 

Fallout from the Fukushima Daiichi Nuclear Power 
Plant (FDNPP) accident in 2011 deposited radioactive 
material over wide areas of North East Japan. Data ob¬ 
tained from air dose rate surveys undertaken after the ac¬ 
cident were primary inputs into the decision to evacuate 
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highly contaminated regions in the early weeks after the 
accident [7] . The data were also used in assessments of ra¬ 
diation doses received by the exposed populations [8, 9]. 

The air dose rate surveys employed various practices 
to gather measurements. In situ monitoring and vehicle 
surveys were conducted at ground level [1, 10, 11]. Heli¬ 
copters and fixed-wing aircraft were used to gather data 
covering wide areas from altitude [7, 12]. The comple¬ 
mentary programs have now provided coverage for much 
of Japan. 

One limitation for converting between air dose rates 
and soil radioactivity levels is that the conversion factors 
are typically derived assuming perfectly flat land sur¬ 
faces [4]. Uneven land topographies can alter air dose 
rates. For example, an ascending slope near to a detec¬ 
tion point will result in a higher air dose rate compared 
to if the ground were flat, as, assuming the total source 
activity and horizontal distribution are unchanged, the 
radioactive source is closer to the detector in the former 
case. The converse is true for slopes descending away 
from detectors; the air dose rate will be lower than if 
there were flat ground as the source is further from the 
detector [13]. 

Topographic effects can also cause uncertainty in con¬ 
verting measurements taken at altitude in aircraft radi¬ 
ation surveys to ground-level values. For instance, the 
ground-level air dose rate can be estimated from a high 
altitude measurement by using a conversion factor be¬ 
tween the 1 m ambient dose equivalent rate measured at 
a calibration site to a detector response measured in an 
aircraft at various elevations above that site [12, 14]. As 
the calibration sites are typically located in flat areas to 
allow vehicle access to make the ground-level measure¬ 
ment, it is not clear a priori that the derived conversion 
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factors will also be applicable if the survey goes on to 
cover mountainous areas [7, 12, 15]. 

Previously Allyson investigated the effect of topogra¬ 
phy on unscattered photon fluence rates at various al¬ 
titudes using an analytical method [16]. He considered 
^^^Cs fallout on the soil surface in valleys and on moun¬ 
tain ridges. The largest effects were seen for steep to¬ 
pographic features and for high detector heights. Saito 
carried out Monte Carlo calculations for various different 
land topographies, but with natural radionuclides, and 
found deviations in air kerma rates of up to 50% [17]. 
Schwarz et al. introduced a method for correcting air¬ 
borne gamma measurements for the perturbations in¬ 
duced by large scale topographic features [13]. 

Satoh et al performed Monte Carlo calculations for 
ambient dose equivalent rates from radiocesium fallout 
for smaller topographic features, in particular 5 m slopes 
up to 40° in inclination [18]. Slopes of this size are simi¬ 
lar to engineered slopes and embankments found within 
Fukushima Prefecture, and they found that topography 
could lead to a 20 % enhancement in ground-level dose 
rates. 

The mountainous region lying to the west of the 
Fukushima Daiichi site, known as the Nakadori area of 
Fukushima Prefecture, suffered large deposits of radioac¬ 
tive ^^^Cs and ^^^Cs. This situation contrasts with the 
1986 Chernobyl accident, where the most serious deposits 
of fallout were over relatively flat areas of Ukraine and 
Belarus in comparison. Therefore we ask how large the 
effect of wide area topographical features, such as moun¬ 
tains and valleys, can be on air dose rates. We consider 
the uncertainty induced by topographic effects when con¬ 
verting 1 m air dose rates into soil activity levels. We 
also consider air dose rates at higher elevations above 
the ground, to analyze the topographic uncertainty when 
converting radiation measurements taken within aircraft 
to ground-level values. 


II. METHODS 



FIG. 1. Angled ground geometry. The detector is the blue 
volume. Fallout is indicated as a red layer on the land surface, 
and soil is beige. 



FIG. 2. Land topography near to FDNPP. Marked are typical 
land inclination angles for the area, as measured over 500 m 
to 1 km distances shown by black bars. Map tile Gopyright 
Thunderforest, map data Gopyright OpenStreetMap contrib¬ 
utors (2014). 


We considered the effect of topography on air dose 
rates at heights h [m] above land surface, from 1 cm up to 
300 m. The higher altitudes cover ranges typically sam¬ 
pled in aircraft surveys. Note 1 m above ground level is 
the standard height for quoting ambient dose equivalent 
rates in the event of a nuclear accident [19]. 

Two idealized geometries were simulated using the 
Monte Carlo radiation transport code PHITS [20]. The 
first was angled ground at the intersection of two planar 
surfaces (fig. 1). The geometry consisted of soil and air 
volumes within a 1 km radius bounding cylinder. Flat 
ground changes abruptly to a uniformly angled slope 
along one direction. The height of the air in the sim¬ 
ulation above the flat land was 400 m, and the minimum 
depth of soil at all locations was 1 m. 

Gamma tracks were detected in 10 m length and 
breadth cuboids above the fiat ground adjacent to the 


slope. To generate sufficient statistics, the vertical thick¬ 
ness of the detector volumes varied between 0.5 cm at 
ground level, to 2 m at the highest altitudes. 

The angle of inclination 0 [°] of the slope was sys¬ 
tematically varied between —20° and 20°. This range of 
land gradients is characteristic of the urban and moun¬ 
tainous areas contaminated with radioactive cesium in 
Fukushima Prefecture, as can be seen from the map in 
fig. 2 that covers the area immediately west of FDNPP. 
Positive angles indicate slopes that ascend up from the 
flat ground, and negative angles the converse. In the 
simulations the angle ^ = 0 corresponds to uniformly flat 
terrain. 

The second geometry considered was a conical geom¬ 
etry. In this geometry, slopes ascend (or descend) away 
from a circular area of flat ground, radius 100 m, in the 
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FIG. 3. Conical geometry with two detector volumes - D1 and D2. 


r [m] 




FIG. 4. Simulation method for the flat ground model. Red indicates radiocesium (the source region), blue detector volumes, 
and beige uncontaminated soil. 


center of the simulation (fig. 3). This geometry approx¬ 
imates detector locations above hilltops and valley bot¬ 
toms. 

Photon detection took place in cylinders above the cen¬ 
ter of the fiat area (denoted Dl), and annular cylinders 
above the fiat ground bordering the slope (D2). The ra¬ 
dius of the Dl cylinders was 10 m, and the horizontal 
thickness of the D2 annular cylinders was also 10 m. The 
vertical thickness of the detector volumes varied between 
0.5 cm and 2 m, depending on the altitude of the dose rate 
point. The slope angle was varied between —20° and 20°, 
and other parameters for this geometry were similar to 
the angled ground simulations. 

The results for these topographies were compared 
against a reference case of uniformly fiat ground. In this 
case there is x-y invariance in the geometry and it is 
possible to rescale the source and detector volumes to 
dramatically improve the computational efficiency of the 
simulation [21] (see fig. 4 and ref. [22]). This method also 
offers the benefit that it allows the contribution of differ¬ 
ent areas of the source on the ground to the overall air 
dose rate to be calculated. In particular we calculated 
the contribution from annular rings of source, between 
inner and outer radii of r [m] and r + 1 respectively, in 
order to check the finite size scaling of the results for the 


air dose rates. 

All models consisted of two materials - air and soil. 
The atomic composition of these substances matched 
the values given in ref. [23]. The density of air was 
1.2 • 10“^gcm“^ and soil was 1.6gcm“^. 

The radiation source was photons with either the ^^^Cs 
or ^^^Cs decay energy spectrum. The spectra were drawn 
from NuDat2, which is based on the Evaluated Nuclear 
Structure Data File (ENSDF) [24]. 

The radiation source was a thin layer of contamination 
on the surface of the ground, or at 1 cm depth within the 
soil. The former case equates to a relaxation mass per 
unit area of = 0 [gcm“^]. 

Note the radiocesium activity per unit area of land 
when projected onto a horizontal plane was constant for 
all simulations. This allows fair comparison of the effect 
of topography on dose rates between the different slope 
inclinations. In other words the total radioactivity of 
cesium source was constant for all slope inclinations in 
the angled ground and conical simulations. 

Ambient dose equivalent rates were calculated by mul¬ 
tiplying the flux detected at each height by conversion 
factors for i^*(10) from ICRP 74 [25]. The flux was 
tallied within 0.5 keV wide energy bins and logarithmic 
interpolation was used between the point i7*(10) conver- 
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FIG. 5. The fractional yield of the infinite half-space ambient dose equivalent rate for a radius r source on uniformly flat 
ground. Colored lines are for ^^^Cs, and superimposed dashed lines for ^^^Cs. 


sion factors. 

Note that none of the models studied here include any 
fine detail for the land structure or use, such as buildings, 
tarmac, surface water, vegetation or local variations in 
topography. 

III. RESULTS 
A. Field of View 

The fields of view for airborne gamma measurements 
can extend to hundreds of meters in radius on the 
ground [16, 26], so first we demonstrate that our simula¬ 
tion spaces were sufficiently large to avoid to finite size 
effects in the calculated dose rates. Fig. 5(a) displays 
the fractional air dose rate for limited radii sources com¬ 
pared to the asymptotic limit for an infinite half-space. 
These simulations were of a thin radiocesium layer on the 
surface of uniformly fiat land. 

The higher the detector elevation, the wider the field 
of view of the detector. At 300 m elevation the dose rate 
from a 900 m radius source is 99.8% of the value for a 
1 km radius source. This indicates that there is no signif¬ 
icant contribution to the total dose rate from any source 
located beyond 1 km. 

The field of view narrows when the radiocesium layer 
is within the soil pack (fig. 5(b)). The effect is most 
pronounced for detectors close to the land surface. This 
is because lower elevation detectors have wider fields of 
view with respect to detector height, i.e. larger solid 
angles, hence the photon path length within soil is 
longer [27]. 

For both surface and 1 cm deep soil layers, the fields of 
view coincide for ^^^Cs and ^^^Cs contamination. This is 
because the i34/i37Qg decay photons that make the dom¬ 
inant contribution to the dose rate are emitted within 
the energy range 0.6 to 0.8 MeV, and so have compara¬ 
ble mean free paths in air and soil. For the angled ground 


and conical geometries, we simulated only ^^^Cs sources 
and infer that the results for ^^^Cs will be quantitatively 
similar. 


B. Angled Ground Geometry 

Slopes ascending up from flat ground result in in¬ 
creased dose rates, as the source is closer to the detector 
than for uniformly flat ground (fig. 6, positive 0). Con¬ 
versely, slopes descending away from the detector result 
in lower dose rates, because the radiation has to travel 
further to reach the detector. 

The effect for 1 m high detectors is relatively small 
when the ^^^Cs fallout is on the soil surface (fig. 6(a)). 
The maximum increase in the dose rate for ascending 
slopes up to 20° is 3%. A change in behavior is visible 
for the descending slopes between —5° and —10° inclina¬ 
tions. The dose rates decrease faster relative to the flat 
ground case compared to the increase for positive slope 
inclinations. This occurs because, for sufficiently large 
negative inclinations, the soil shields the direct path be¬ 
tween radiation originating on the slope and part of the 
detector volume. 

For higher altitude detectors the effect of the topog¬ 
raphy on dose rates is generally stronger, particularly at 
300 m elevation where the dose rates change by up to 
26% relative to the flat ground case. However, for the 
higher elevations, the geometries mean that slopes de¬ 
scending away from flat ground do not shield the source 
from the detector. For 50 and 300 m detector elevations 
there is little difference between radiocesium layers on 
the soil surface and at 1cm deep (c.f. 6(a) and (b)). 

The effect of topography for 1 m detectors is slightly 
stronger for ascending slopes when the contamination 
layer is within the soil than when it is on the surface. 
Conversely for descending slopes, there is a smaller effect 
of topography as the narrower field of view means that 
attenuation by the slope itself is of lesser importance. 
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FIG. 6. Relative change in ambient dose equivalent rates for different land inclination angles in the angled ground geometry. 
Schematics of the geometries, sources and detector locations for positive/negative 0 are shown below the abscissa. 



FIG. 7. Ambient dose equivalent rate conversion factors as a 
function of elevation above perfectly flat ground. 


A limitation of these results is that they only apply 
for the simulated detector locations. It is possible to 
determine qualitatively how the dose rates would change 
at different positions. If the detector moves away from 
the slope and towards the flat region, the limiting dose 
rate is identical to the dose rate for uniformly flat ground. 

Moving in the opposite direction, i.e. up or down the 
slope, the limiting value is the dose rate at height h cos 0 
above a uniformly flat surface. Thus there is a change in 
the effective height of the detector. This is a purely geo¬ 
metrical effect attributable to the difference between the 
vertical height above the slope and the shortest distance 
to it. 

To elucidate the consequence of this effect, fig. 7 shows 
i7*(10) conversion factors per unit ^^^Cs and ^^^Cs sur¬ 
face activity for various elevations above flat ground. The 
ambient dose equivalent rate shows a logarithmic reduc¬ 
tion with height. Therefore, the effective reduction of the 
detector altitude by cos 0 results only in a small increase 
in the air dose rate. 


C. Conical Geometry 

The second type of land geometries considered were 
conical slopes surrounding flat land. The inclination an¬ 
gle of the slopes was systematically varied, and the air 
doses rates compared against the uniformly flat ground 
reference case, fig. 8. 

The trend of the data is similar to that for the angled 
ground geometry, however the magnitude of the effect is 
stronger as slopes now surround the detectors. Slopes as¬ 
cending near to a detector result in higher air dose rates, 
and slopes descending away lead to smaller dose rates. 
The effect is most pronounced for the largest slope incli¬ 
nations, and for detector locations D2 that are adjacent 
to the slope. The dose rates differ by more than 35 % for 
detectors at 300 m elevation and large inclination angles 
when compared to perfectly flat land. 

There is comparatively little difference between the two 
different source locations for 50 and 300 m high detectors, 
c.f. figs. 8(a) and (b). For the Im ambient dose equiva¬ 
lent rates, there are slight variations between the results 
for the surface and 1 cm deep source layers. These vari¬ 
ations are attributable to the different fields of view and 
consequentially the relative effect of bringing the source 
closer to or further from the detector by changing the 
slope inclination. However, the magnitude of the effect 
is always small, and the dose rate is always within 23 % 
of the dose rate for perfectly flat land. 


IV. DISCUSSION 

A. Ground-Level Air Dose Rates 

The simulations were based on large topographic fea¬ 
tures, such as mountains and valleys, and the studied 
land inclination angles varied between —20° and 20°. 
There are only small deviations in the 1 m ambient dose 
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FIG. 8. Relative change in ambient dose equivalent rate with slope inclination for the conical geometry. 


equivalent rate compared to the fiat ground case for this 
inclination angle range. Larger deviations are possible 
for smaller but steeper topographic features such as en¬ 
gineered slopes [18]. 

The difference between the air dose rates for uneven 
and perfectly flat ground sets the uncertainty due to to¬ 
pography when converting 1 m dose rates into soil activ¬ 
ity levels using standard conversion factors [4] . If the air 
dose rate is higher than otherwise would be the case for 
flat land due to the presence of ascending slopes near to 
the detector, the estimate for radiocesium concentration 
within the soil will be an overestimate. The soil activity 
will correspondingly be underestimated if slopes descend 
away from a detector location. 

These results apply under the assumption of spatially 
homogeneous radiocesium fallout. In practice there can 
be large spatially variability in soil activity levels over 
relatively small areas [26] . As the field of view for 1 m 
ambient dose equivalent rates is on the order of 10 to 
100 m, uncertainty in predictions for soil activity is more 
likely to be due to spatial heterogeneity of contamination 
than topographic effects. 


B. Airborne Dose Rate Surveys 

The effect of topography on air radiation dose rates 
was generally stronger for higher altitude detector loca¬ 
tions, particularly at 300 m, than for detectors close to 
the land surface. If converting between aerial gamma 
measurements and ground-level dose rates using a con¬ 
stant factor derived for flat land [14], the uncertainty for 
uneven topographies may be as large as 60 %. However, 
as this represents the largest deviation seen in the simu¬ 
lation results, the topographic uncertainty is likely to be 
smaller than this value in practice. 

The effect of the topographic uncertainty on the 
ground-level predictions from airborne surveys employing 
conversion factors suited for flat land is as follows. Pre¬ 


dictions for ground-level dose rates near hilltops will be 
underestimates, and correspondingly predictions made 
for near valley bottoms will be overestimates [16]. 

V. CONCLUSIONS 

For the geometries and detector locations studied here, 
the air dose rates for varying land topographies are al¬ 
ways within 60 % of their values for perfectly flat land. At 
detector heights near to ground level, the effect of large 
topographic features on air dose rates was small. Other 
factors, such as spatial heterogeneity of the contamina¬ 
tion, land use, buildings and decontamination efforts, will 
be more important sources of uncertainty in calculations 
for soil activity levels based on measured dose rates. 

The largest effects of topographic uncertainty were 
seen for the high altitude detectors. Wide-area topo¬ 
graphic features, such as mountains or valleys, could 
therefore be significant sources of uncertainty when in¬ 
terpreting the results from airborne dose rate surveys. 
Mountainous areas are also the regions where airborne 
gamma surveys are most powerful, as it is more difficult 
to access mountainous terrain at ground level to perform 
radiation surveys. 

It is worth noting that comparable levels of topo¬ 
graphic uncertainty were found for 50 m elevation and 
ground-level detectors. Unmanned aerial vehicles can op¬ 
erate at lower altitudes than traditional fixed-wing air¬ 
craft or helicopters [15, 28], and therefore drone surveys 
will be less susceptible to topographic uncertainties in 
their processed output. 
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